Precision measurement of spin-dependent interaction strengths for spin-1 and spin-2 8/ Rb atoms 
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We report on precision measurements of spin-dependent interaction-strengths in the 87 Rb spin-1 and spin-2 
hyperfine ground states. Our method is based on the recent observation of coherence in the collisionally driven 
spin-dynamics of ultracold atom pairs trapped in optical lattices. Analysis of the Rabi-type oscillations be- 
tween two spin states of an atom pair allows a direct determination of the coupling parameters in the interaction 
hamiltonian. We deduce differences in scattering lengths from our data that can directly be compared to theo- 
retical predictions in order to test interatomic potentials. Our measurements agree with the predictions within 
20%. The knowledge of these coupling parameters allows one to determine the nature of the magnetic ground 
state. Our data imply a ferromagnetic ground state for 87 Rb in the / = 1 manifold, in agreement with earlier 
experiments performed without the optical lattice. For 87 Rb in the / = 2 manifold the data points towards an 
antiferromagnetic ground state, however our error bars do not exclude a possible cyclic phase. 

PACS numbers: 03.75.Lm, 03.75.Gg, 03.75.Mn, 34.50.-s 



Contents 



Ctf 



T3 

o 
o 



> 

00 

m 

o 

o 

^D 
O 
-i— > 



i 

-a 

c 
o 



X 



I I. Introduction! 

I II. The ory of coherent collisional spin dynamics! 



^^^rjin^hangin^o llisions in optical lattices! 



I^JwoJeve^^sterrjL 



C. Three level system! 



IlII. Experimental sequence! 
I IV. Experimental results! 



^^Dvjiamk^rijhe^^^ljiy^erfine^tate 



B^^vmmik^njhe^^^^rvjjerfme^tate 



C 1 ^e^ermjrrmg_s£atterin g length differences! 



IX Damping mechanisms! 
I V. Nature of the magnetic ground state! 
IVI. Conclusion! 
I Acknowledgments! 

I Coupling strengths! 

I Rabi dynamics! 

I References! 



I. INTRODUCTION 



9 
9 

10 



The experimental realization of Bose-Einstein condensa- 
tion (BEC) in optical traps 10, H H H has opened the 
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possibility to investigate novel effects not accessible in mag- 
netic potentials, where only specific atomic spin states can 
be trapped. In an optical trap the spin-degree of freedom is 
liberated, allowing to study intriguing phenomena in multi- 
component gases such as spin dynamics OHyll> spin waves 
O lal, spin-mixing I9l llOll and spin-segregation [11]. The 
mechanism behind these phenomena is a coherent interac- 
tion between two atoms which changes the individual spins 
of the colliding particles while preserving the total magneti- 
zation o eei mm. 

In the low-energy limit, interatomic interactions are char- 
acterized by the corresponding interaction strength, usually 
parametrized by the s-wave scattering length a 1171 Hal . In 
general, this scattering length depends on the spin states in the 
incoming and outcoming scattering channels. For two collid- 
ing spin-/ particles, the interaction strength is characterized 
by 2/ independent scattering lengths Of, labelled by the total 
coupled spin F of the two particles. For bosons (fermions) 
the total spin F takes even (odd) values between and If, 
This spin-dependence of the interaction potential is reflected 
in a many-particle system by a spin-dependent interaction en- 
ergy. The spin configuration which minimizes this energy at 
zero magnetic field is energetically favoured and is called the 
magnetic ground state of the system. 

For interacting spin-1 gases this magnetic ground state 
can be either ferromagnetic or antiferromagnetic in nature, 
i.e. spin dependent interactions favour the relative orienta- 
tion of two atomic spins to be either parallel or antiparal- 
lel I12l Il3l . In the spin-2 case an additional cyclic phase 
can also arise 1141 Il5l ll9B . Here spin-spin correlations are 
revealed for three particles, whose spin orientations tend to 
form an equilateral triangle in a plane. These systems are 
prom ising in order to study quantum magnetism phenomena 
\Ya IT3L !l4l Hal , and to create strongly correlated magnetic 
quantum states in optical t raps 1151 l2CLl2lLl22i 12311 and optical 
lattices JMillllllllMlMlMlIinillll. 

Many experiments on multi-component quantum gases in- 
volve 87 Rbl3ll4ll5l l34ll . the "workhorse" of present day BEC 



experiments. Several experiments have investigated the dy- 
namics of a spin-1 87 Rb Bose-Einstein condensate (BEC) in 
order to identify the magnetic ground state |J, |j] . There, a 
slow spin dynamics has been observed on a second timescale 
leading to a final state close to the predicted ferromagnetic 
ground state at zero external magnetic field. Recently, an ex- 
periment observing coherent spin dynamics in a BEC 1 34113511 
has shown the spin-1 magnetic ground state of 87 Rb to be 
ferromagnetic. For 87 Rb spin-2 atoms the magnetic ground 
state has been predicted to be antiferromagnetic, but close to 
the phase boundary of the cyclic phase. Recent experiments 
1 3, 5], again observing the spin dynamics for various popu- 
lation of the Zeeman sublevels, seem to agree with this pre- 
diction. One possible difficulty, however, arises from the fact 
that even small magnetic fields can perturb the system and 
pin it to some configuration which is not the ground state at 
zero external field. Another difficulty in precisely determin- 
ing the ground state stems from the fact that the differences of 
the scattering lengths driving the spin dynamics is very small 
(typically only a few percent of the bare scattering lengths), so 
that the time scale on which the system relaxes to its ground 
state is very long. More recently, a novel method has been 
proposed theoretically to identify the magnetic ground state 
of 87 Rb spin-2 atoms 1 36] circumventing many of the prob- 
lems described above by observing the time evolution of an 
initial state with specific population and phases of the Zeeman 
sublevels. 

In this work, we choose another approach and present a 
method based on the coherence of the spin evolution to di- 
rectly deduce the absolute value and sign of the various spin- 
dependent interaction terms in the hamiltonian. From the 
measured oscillation frequencies we are able to infer the coef- 
ficients in the interaction hamiltonian and value of the differ- 
ences in scattering lengths. The knowledge of these allows in 
turn to determine the magnetic ground state. There are several 
advantages of this method. First, our system is intrinsically 
free from mean-field shifts. Second, the optical lattice leads to 
a significant enhancement of the coupling strength determin- 
ing the time scale of the spin-changing dynamics. Therefore, 
even for the spin-2 case, where strong losses are expected, co- 
herent spin-changing collisions can be observed. This allows 
for a direct comparison with calculations based on a recent 
prediction for 87 Rb ground state s-wave scattering lengths. In 
this sense our measurement constitutes a stringent test of the- 
oretical interatomic interaction potentials. 

This paper is structured as follows. Section UTI briefly re- 
calls the fundamental mechanism of spin changing collisions 
and introduces the collisionally induced Rabi-oscillations be- 
tween two-particle Zeeman states in an optical lattice. We 
give a detailed description of our experimental sequence in 
Section [H]] In Section HV1 we present the main experimental 
results and deduce the relevant parameters characterizing the 
spin-dependent interaction. In Section|V]we describe how the 
precise determination of the collisional coupling constants can 
identify the magnetic ground states of spin-1 and spin-2 87 Rb 
ground state atoms. 



II. THEORY OF COHERENT COLLISIONAL SPIN 
DYNAMICS 

A. Spin changing collisions in optical lattices 

Let us consider two colliding spin-/ 87 Rb atoms initially 
in the single-particle Zeeman states with magnetic quantum 
numbers m\ and W2. The particles can undergo a spin chang- 
ing collision and be tranferred into a final state character- 
ized by the quantum numbers mj and m\. In binary colli- 
sions between alkali atoms the total spin pr ojection on the 
quantization axis is conserved I^ fl0l[l2l[l3ll . which implies 
m\ + mi - m-i + m\. This limits the number of accessible final 
states to those that have the same total magnetization as the 
initial one. The interaction can be described by a potential of 
the form 1 12] 



2/ 



V{A - A) = v s 6(?i - n) = I^sf'Pf 5{A - n) 



(D 



F=0 



where gf - (47tH 2 /M) ap, M is the mass of one atom, Pf is 
the projection operator onto the subspace of total spin F, ap 
is the s- wave scattering length for two atoms colliding in such 
a channel with total spin F, and Ti - h/2n, with h Planck's 
constant. For bosons, the total spin F is restricted to even 
values only 1 12, 19]. 

In a deep three-dimensional periodic potential (see Section 
nil , the situation can be adjusted such that at many (~ 10 5 ) 
potential minima (lattice sites) there are two atoms trapped in 
the vibrational ground state of the potential. Throughout the 
following discussion we assume a sufficiently deep potential 
such that the atoms can neither leave the trap by tunneling nor 
be excited via interactions into a higher vibrational state. Con- 
sequently, the atoms remain in the lowest vibrational state at 
all times. Because the spatial dynamics is frozen, the ensuing 
dynamics caused by spin-dependent interactions only affects 
the spin component of the two-particle wavefunction. To de- 
scribe this collisional dynamics, we write two-particle atomic 
wavefunctions in the form 



W) = \fo)\\hh* S\f, m;f,m'). 



(2) 



Here 1,2 labels the two atoms in the pair, S is the symmetriza- 
tion operator, |0o) is the spatial wavefunction corresponding 
to the vibrational ground state in each well. Also \f,m;f,m') 
denotes the (unsymmetrized) two-particle spin wavefunction. 
In the following we describe the atom pair by its spin wave- 
function only and imply a constant spatial wave function |0o) 
that is not explicitly denoted. With this notation, the inter- 
action potential Q connecting an initial two-particle state 
|/, «i;/, TO2) with a final two-particle state \f,my,f,rri4) can 
be written as 



(f, my, f, m 4 \V\f, m, ; /, m 2 ) = Ux Aa^Z 



(3) 



where the two quantities on the right hand side are given by 

A<£ = Zlio ZL- F (f> m; f, m 4 \F m) a F (4) 

U =±f x jd 3 r \<p \ 4 - (5) 



Here, Aa^lmt is given by a difference of the scattering lengths 
a F weighted by the relevant Clebsch-Gordan coefficients 
(f,my,f,m4\Fm) connecting the initial and final spin states 
through an intermediate, coupled spin state \Fm). The pre- 
factor U contains all constants and the spatial wavefunction 
overlap. It can be calculated using 



The effective Rabi frequency Q.' if is given by 



U - Uno + f/anharm + C/fe 



(6) 



Uho is calculated by approximating the lattice potential by 
a harmonic well near each minimum, and by using the low- 
est harmonic oscillator eigenfunction as <po Km . In addition, 
we include two corrections. First, using Wannier functions 
from a band structure calculation for the wave function <po in- 
stead of a gaussian yields a negative correction term f/anharm 
which is on the order of 5%. Second, at short distances the 
spatial wave function of the colliding particles is expected to 
deviate from the product form in (0- A preliminary theo- 
retical investigation |38] indicates that this deviation may in- 
troduce a negative correction Ui nt . We have experimentally 
performed a consistency check for the prefactor U by mon- 
itoring the collapse and revival of the interference pattern as 
performed in 1 39]. We measure the interaction matrix element 
U for the spin state |/ = 1, m = -1), Assuming a known scat- 
tering length for this spin state we infer an approximate value 
of U consistent within 10% with calculated values including 
only the correction due to the anharmonicity Oanharm- Hence, 
the actual values we use in the following include the first cor- 
rection term f/ a nharm, but not the second U[ nt . The values are 
U = 2n x (30.4 ±1.1) Hz/a B for a lattice depth of 40 E r , and 
U = 2tt x (33.0+ 1.4) Hz/aB for 45 E r . Here a B is the Bohr ra- 
dius, E r is the single photon recoil energy E r = h 2 /2MA 2 , with 
A the wavelength of the lattice laser. The error is given by our 
uncertainty in optical potential depth, for which we assumed 
an upper bound of 10%, but a possible additional systematic 
error of U due to the correction {7; nt is not explicitely given. 
Finally it should be noted that the collapse and revival experi- 
ment has been performed at a magnetic field of B - 150 G in 
the magnetic trap, whereas the dynamics driven by the spin- 
dependent interaction occurs at low magnetic fields around or 
below 1 G. 



B. Two level system 

In the experiments described below (see cases a-c in Fig- 
ure [Q, we prepare all atom pairs in the lattice in a specific 
spin state |^,-). This initial state is in general coupled to one 
or several final states |<A/) by the spin-changing interaction. 
In the simplest cases, where the atom pair is initially pre- 
pared in |1, 0; 1, 0> (Figure[TJi) or in |2, -1; 2, -1) (Figured), 
the system can access only one final state, 1 1,-1-1; 1,-1) or 
|2, 0; 2, -2), respectively, while conserving total magnetiza- 
tion. 

As for any two-level system, the probability to find the atom 
pair in the final state is given by the Rabi formula 



Q 2 

2Q' 2 . 



[l-cos(Q.' if t)]. 



nJ/=V^+4- 



(8) 



where the coupling strength Q,y is determined by the ma- 
trix element (if/f\Vs\tffi) of the spin-changing interaction, and 
where the detuning Sif can be decomposed into two parts: 



8 if - £int + S B . 



(9) 



The first detuning 6- mt arises from the difference of interaction 
energies in the initial and final states (iffi\Vs\iffi) - (if/f\Vs\if/f). 
The second detuning 6 b is due to the energy shift of initial and 
final states in a non-zero magnetic field. Since the total mag- 
netization is conserved, the system is not sensitive to first or- 
der Zeeman shifts, because the initial and the final states expe- 
rience the same first order Zeeman effect |40]. Consequently 
the detuning is to leading order determined by the quadratic 
Zeeman shift. Therefore the term 6b in Equation (|5Jl is given 



by the second-order Zeeman-effect 6b • 



(MbB) 1 



where Whf S de- 



(7) 



notes the hyperfine splitting. The corresponding detuning <Sb 
is shown versus magnetic field in Fig.^1. 

As indicated in the previous section, the two parameters of 
the Rabi model can be written in the form Aa x U, where the 
value of Aa, different for Q,y and Sif, can be calculated for 
each specific case (see Table IIHI . It should be noted that for 
87 Rb the bare scattering lengths af are almost equal; recently 
calculated values |41] are ao = (101.78 + 0.2) ob and ai — 
(100.40 + 0.1) a B for spin-1. Predicted values for spin-2 O 
are a Q = (87.93 + 0.2)a B , a 2 = (91.28 + 0.2) a B and a A = 
(98.98 ± 0.2) a B , where the Bohr-radius a B = 5.29 x 10" 11 m. 
Thus the differences in scattering lengths describing the spin- 
dependent coupling give values of Aa on the order of a few 
percent of the bare scattering lengths af. 

Since the dependence of the magnetic field dependent de- 
tuning with external field is known, a measurement of the ef- 
fective Rabi frequency Q,' if for different magnetic fields allows 
to extract the absolute value of the coupling strength O,/ and 
both the absolute value of the constant detuning 6 mt and its 
sign. The form of Equation l|S} implies that if 5; nt and 6b have 
the same sign, even at zero external magnetic field Rabi os- 
cillations do not show unity contrast, because the constant 
detuning 6 mi is present at all times. This can be overcome 
e.g. applying a microwave field which can shift the two two- 
level states with respect to each other in energy. The effect of 
this AC Zeeman effect on spin dynamics in an optical lattice 
has been described in detail in |42]. 



C. Three level system 

The situations investigated in the previous section were de- 
scribable as pure two-level systems. Another possibility (case 
(c) in Figure ^ is mat t\v° final states \4>j\) and \4>f 2 ) fulfill 
the condition of a conserved total magnetization. The dynam- 
ics can then be modelled by a coupled three-level system, as 
has been calculated in detail in |43]. Assuming that the direct 



The eigenfrequencies then read 



f=2 

m f -2 

f=1 



-1 +1 +2 
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i(< 



6\ + fo + Q cos - 1 



£ 

3 

l / - In-l 

CO - -l(5i +fo + Q cos — - — 

1 / - 2ti + £ 

co - — \6\ + 62 + O cos — - — 



(12) 
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f=2 
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where 
Q 



In - arccos t=- [9(Oj + Qj) 

-4(25, - <5 2 )(2<5 2 - *i)](*i + £2) - 27Q^ 2 ] . (13) 



0.5 1.0 1.5 2.0 
Magnetic field B (G) 



In the case where |</</,) is separated by a very large energy dif- 
ference from the other two states, 62 » Si, the system effec- 
tively reduces to a two-level system oscillating between |(//,-) 
and |i///,). The relevant oscillation frequency then becomes 

lim ( 5,^ oo ((y - tiT) - J& 2 + 6 2 . In this limit, a description as 
a two-level system \\jfi) <-> |^«) is valid. 

Comparing to the situation depicted in FigureQJ this means 
that for low magnetic fields (B < 0.6 G) a model including all 
three states has to be used in order to describe the dynam- 
ics. For increasing magnetic fields, the second order Zeeman 
shift detunes the final state |0y,) = |2, +2; 2, -2) faster than 
the other final state \ij/f x ) - |2, +1; 2,-1) (see Figure^!), such 
that for large enough magnetic fields (B > 0.6 G) a two-level 
model only including two states |2, 0; 2, 0) and |2, +1; 2, -1} is 
sufficient. 



FIG. 1 : (a-c) Shift of energy levels due to second order Zeeman ef- 
fect in the 87 Rb ground state with spin-1 and spin-2. The various 
processes showing the spin changing dynamics between two-particle 
Zeeman states are indicated in different colors, (d) Corresponding 
magnetic field dependent detuning Sb for selected combinations of 
two-particle states versus magnetic field. 



coupling between \t//j) and \\j/^) is negligible, the correspond- 
ing interaction hamiltonian reads 



h 

V= 2 



( Qi 
Qi 28] 
Q.2 




Q 2 

26 2 ) 



(10) 



where Q, and 6j describe the coupling strength and total detun- 
ing, respectively, and i = 1(2) labels the process |t/r,) «-> |^/,) 
(|^/j) «-» \i/ff 2 )). Similar to the treatment in |44] we find the 
eigenfrequencies co as solutions to the secular equation 



to 3 - co 2 {5\ + 62) + co I 6162 



<$ 



- Q i 



"I 



S 2 



0. (11) 



III. EXPERIMENTAL SEQUENCE 

In order to realize the different situations of spin-dynamics 
in an optical lattice as depicted in Figure [Oa-c), we start 
our experimental sequence by creating a 87 Rb Bose-Einstein 
Condensate (BEC) containing around 2 x 10 5 atoms in the 
\f = \\nif - -1) hyperfine state in an almost isotropic mag- 
netic trap (trap frequency com ~ 2nY. 15 Hz). We then load the 
sample into a combined magnetic and optical trap by slowly 
increasing the power of three mutually orthogonal standing 
waves to a final trap depth of either 40 E r or 45 E r , driving 
the system into the Mo tt-insulat or regime [ 58] . Here E r is the 
single photon recoil energy E r - h 2 /2MA 2 , with A - 840 nm 
the lattice laser wavelength. In the later discussion of system- 
atic errors we assume an uncertainty in the potential depth of 
~ 10%, which should be understood as an upper bound. The 
particular form of the lattice-intensity ramp has been taken 
from 1 45], minimizing excitations in the system |46]. It has a 
total duration of 160 ms (see Fig.|2l. 

In order to observe spin oscillations, the magnetic trap is 
switched off, leaving the sample in a pure optical trap. The 
spin-polarization of the atoms is preserved by a small homo- 
geneous offset-field of ~ 1 .2 G which is switched on 5 ms 



A Magnetic fields 




f=1 




f=2 




gion around the corresponding atom cloud (see dashed boxes 
in Figures |2),c). The total atom number A^t is counted in 
a large region including all magnetic sublevels. However, the 
different atom clouds corresponding to different magnetic sub- 
levels are not well enough separated. In order to avoid dou- 
ble counting of atoms, we chose counting regions that do not 
overlap, but slightly underestimate the population N m of each 
cloud. Moreover, an additional background signal on the or- 
der of a few percent is counted in the global region used to find 
the total atom number iV tot . The populations in different Zee- 
man states as given later therefore sum up to a value slightly 
(< 5%) smaller than 1. 

Due to the gaussian shape of the lattice laser beams, the 
atoms in the optical lattice experience an additional smooth 
confinement and therefore the system favours the formation 
of Mott shells with distinct, integer number of atoms per site 
I47tl48ll49ll . For our trap parameters we calculate from a mean 
field model that a central core of atom pairs contains approx- 
imately half of the atoms 1 50 ] . This core is surrounded by a 
shell with isolated atoms, whereas the number of sites with 
more than two atoms is negligible. Due to this shell structure, 
not all atoms contribute to the spin oscillations. Consequently 
the measured amplitude of the spin oscillation - normalized to 
the signal from all atoms - is decreased compared to the value 
predicted by the Rabi model. The measured relative amplitude 
is expected to be 



N +1 + N- 
Mot 



= nP f 



(14) 



where n is the fraction of atoms in sites with double filling. 



FIG. 2: (a) Sketch of the experimental sequence, (b) Typical camera 
image for atoms in / = 1 and (c) / = 2. The dashed boxes illus- 
trate areas in which the atom number for each magnetic sublevel is 
counted. 



before switching off the magnetic trap (see Fig0. After a 
hold time of 60 ms during which the bias field of the magnetic 
trap decays to below 1 mG, the atoms are prepared in their 
initial state for spin dynamics by application of microwave 
pulses with a frequency around 6.8 GHz and a pulse duration 
between 6 yt/s and 150//S, depending on the specific transition. 
Directly after the last pulse, the homogeneous offset field is 
ramped to its final value. The time constant of the magnetic 
field change is on the order of a few hundred //s; this magnetic 
field has been calibrated by microwave spectroscopy and is 
known to a value better than ~ 10 mG, which in the follow- 
ing we consider as an upper bound of the systematical error 
in the magnetic field. After a variable hold time t the optical 
trapping potential is switched off and the atoms can expand 
during 7 ms of time-of-flight (TOF), while the homogeneous 
magnetic field is still present. During the first 3 ms of TOF 
a gradient field is switched on in order to spatially separate 
the different magnetic sublevels. The population N m of each 
magnetic sublevel m is recorded after TOF with standard ab- 
sorption imaging by counting the fractions N m in a small re- 



IV. EXPERIMENTAL RESULTS 

We have experimentally investigated the spin dynamics in 
the three situations depicted in Figure Qa-c), where the sys- 
tem is initialized in different internal atomic states. The 
initial states for those cases are (i) |/=l,m/ = 0), (ii) 
1/ = 2, nif - -1), and (iii) \f = 2, nif - 0). In most cases, the 
assumption that for an initial state only one final state is ac- 
cessible after a spin changing collision is valid. In the lower 
/ = 1 hyperfine manifold, this is naturally fulfilled, because 
only one combination of two-particle states conserves total 
magnetization. For the upper / = 2 hyperfine state, there are 
in general more initial and final states accessible. For case (ii) 
the description as a two-level system is valid at all parameters 
used in our experiments. For case (iii), as already discussed in 
Section HTCl a sufficiently large external magnetic field can be 
used to tune the system via the quadratic Zeeman-shift into a 
regime where the approximation by a two-level system is suf- 
ficient. For small magnetic fields, however, a description as a 
three-level system must be considered. 

Another difference between the / = 1 and / = 2 cases is 
the role of losses. Due to our preparation in the optical lat- 
tice, three-body recombination processes are suppressed, be- 
cause most atoms are either in singly or doubly occupied lat- 
tice sites. In the / = 1 state, two-body loss rates are very 
low, so that atom losses are negligible. This is not the case 



in / = 2 where hyperfine changing collisions have been seen 
to be significant Onll 15211 and to lead to a dephasing of the 
spin oscillations |:53|. We still observe several coherent spin- 
oscillations in the upper hyperfine state, because of the rela- 
tively large spin-dependent coupling strengths in / = 2. 

In each case, the frequency of the coherent spin oscillations 
can be changed through the magnetic field dependent detun- 
ing 5b- We record the population in each Zeeman sublevel as a 
function of time for the three initial states introduced above for 
different values of the external magnetic field and extract the 
corresponding oscillation frequency. For each process we fit 
the data to the expected behavior of the oscillation frequency 
(0. Since the spin-dependent dynamics is well described by 
one (two) differences in scattering lengths for the / = 1 (/ = 
2) case, we fit the oscillation frequency vs magnetic field be- 
haviour with only one (two) free parameters corresponding to 
the scattering length differences. In order to do this, we make 
use of the relations given in Table [HI] re-parametrizing the 
coupling strengths Q,y in terms of the detunings 6\ nl ; hence, 
we use the fitting parameter 5\ for / = 1, being the interac- 
tion detuning <5j nt for the process |1,0;1,0) <-> |1,+1;1,-1). 
For / = 2 we use the two interaction detunings 62,1 for the 
process |2, 0; 2, 0) <-> |2,+1;2,-1) <-> |2, +2; 2, -2) and 62,2 
for |2, — 1; 2, — 1) <-> |2, 0; 2, -2) as fitting parameters. We do 
this instead of directly fitting the scattering length differences, 
because the factor U linking the resulting detunings £j nt to the 
scattering length differences Aa is not error-free and might be 
systematically shifted. Finally, from the extracted detunings 
6 mt and the calculated values for U (see Section llTA"l . we give 
values for the scattering length differences calculated from the 
extracted detunings <5; nt . 
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A. Dynamics in the / = 1 hyperfine state 



FIG. 3: (a) Spin population oscillations in / = 1 between |1,0; 1,0) 
(O) and |1,+1; 1,-1) (•) at a magnetic field of 0.28G and a lattice 
depth of 45 E r . A fit to a damped sine (solid lines) yields an oscilla- 
tion frequency of 40.2(3) Hz and a damping rate of 7 = 3(2) Hz. (b) 
Measured oscillation frequency of the |1,0;1,0)«-»|1,+1;1,-1) col- 
lision process versus external magnetic field B. The solid line is a fit 
to the expected behaviour according to the effective Rabi-frequency 
il' if . The dashed line is the calculated curve based on the predicted 
scattering lengths |41]. The shaded region reflects our uncertainty 
in magnetic field of ±10mG and the uncertainty in potential depth 
(upper bound of 10%). 



The collisional coupling between two / = 1 atoms is an 
almost ideal example of the Rabi model introduced above. 
The only possible process conserving total magnetization is 
1 1,0; 1,0) <-> |1,+1; 1,-1). In order to observe spin dynam- 
ics, the atomic sample is prepared by a first microwave pulse 
in the \f - 2; mj - 0)-state and immediately a second pulse is 
applied, transferring the atoms into \f - \,mj - 0). A typi- 
cal oscillation for a lattice depth of 45 E r at a magnetic field 
of B — 0.28 G is shown in Fig. [3^. Spin oscillations have 
been observed for various magnetic fields between approxi- 
mately 170 mG and 600 mG For larger magnetic fields, the 
amplitude is strongly suppressed, because the magnetic field 
dependent detuning 5b becomes large compared to the cou- 
pling strength. In this limit, spin oscillations become difficult 
to resolve. The extracted interaction detuning (see Section lvH 
is dimif = 1) = Si = 2nx (-11.8±°° 4 )Hz, where the upper 
error gives the statistical uncertainty, and the lower error the 
systematic uncertainty in lattice depth. 



B. Dynamics in the / = 2 hyperfine state 

For the spin-2 case, coherent spin oscillations driven by in- 
teratomic collisions in an optical lattice have been presented 
in 1 53]. Here we re-analyse the data, including a description 
by a coupled three-level model. 

a. Two-level system. For the case where the dynamics 
starts from |2, — 1;2, — 1), the only possible interconversion 
process that preserves total magnetization is |2, — 1; 2, -1) «-> 
|2, 0; 2, -2). Thus this situation can also be described as a 
two-level system. 

b. Three-level system. For the case where the system is 
initialized in |2, 0; 2, 0), there are, however, two final states 
|2, +1;2, -1) and |2, +2; 2, -2) accessible. The direct pro- 
cess |2, 0; 2, 0) <-> |2, +2; 2, -2) is strongly suppressed, but 
this final state can be populated through a two-step process 
|2,0;2,0) <-> |2,+1;2,-1) <-> |2,+2;2,-2). For sufficiently 



large magnetic field, the state [2, +2; 2, -2) is far detuned so 
that the situation is effectively described by a two-level system 

10. 

For lower magnetic fields, a three-level description as in 
Section lncl must be used in order to describe the population 
oscillation. In Figure^ we show a population oscillation for 
a lattice depth of 40 E r and a magnetic field of B — 0.24 G 
which clearly differs from the simple sinusoidal case. Here, 
the solid lines are a calculation from a three-level model which 
shows good agreement with the measured data. For this cal- 
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FIG. 4: (a) Measured population dynamics for the coupled three- 
level system |2,0;2,0) <-> |2,+1;2,-1) «-» |2, +2; 2, -2) at a mag- 
netic field of B = 0.24 G and a lattice depth of 40 E r . One clearly 
sees the deviation from the pure sinusoidal two-level oscillations and 
the rather large population of |2, +2; 2, -2). The solid lines are cal- 
culations based on a coupled three-level system ilOi with parameters 
corresponding to the experimental realization, (b) Oscillation fre- 
quency vs magnetic field for spin dynamics starting from |2, 0; 2, 0). 
The red line is a fit to the analytical solution of the three-level model 
I12L the dashed line shows the result of a two-level model as pre- 
sented in 1 53]. 

culation we use the fitted values of the two detunings 62,1 and 
62,2 (see below) together with the form of the Rabi parameters 
given in Table |III] With this we calculate the values of the 
three-level model {lOj D. { - 2n x 265,Hz, 6\ - 2k x 33 Hz 
for the coupling of |2,0;2,0> with |2, +1;2,-1), and Q. 2 = 
27r x 134 Hz, 62 — 2jt x 100 Hz for the coupling between 
|2, +1;2, -1) and |2, +1,2, -2). The phenomenological damp- 



ing rates for both processes have been set to y — 30 s -1 as 
measured in |53] and the decay from |2, +2; 2, -2) has been 
neglected. The amplitude has been set to account for ~ 50% 
of the atoms contributing to the spin-dynamics, and the offset 
for each population (see Section lTV D> has been set in order to 
match the data. 

Although the coupled three-level system defined by d!0> 
can in principle be solved, we were not able to fit its solution 
reliably to the data because of the large number of parame- 
ters involved. Instead, we fit the mj - population with the 
same sinusoidal form as used in the previous cases to approxi- 
mate the predominant oscillation frequency. We have checked 
that a Fourier analysis yields oscillation frequencies in agree- 
ment with those found by this fitting procedure. The observed 
oscillation frequency versus magnetic field is shown in Fig- 
ure 0J>. In order to include the three-level character, we fit 
the data for this process with the analytical solution of the 
three-level model. The fit is shown as solid line in Figure 
0Ji, including the dependency of this process on the other pro- 
cess |2, — 1; 2, -1) <-> |2, 0; 2, -2) via the interaction detuning 
#2,2. The theoretical prediction from a two-level system (as 
presented in 1 53]) is shown as dashed-dotted line and does 
not describe the dynamics well for small magnetic fields. A 
theoretical prediction including the third level, however, does 
reproduce the correct behaviour (dashed line in Figure0J)). 



C. Determining scattering length differences 

For the analysis it is important to notice that the dynamics in 
a spin dependent collision of two spin-1 (respectively spin-2) 
atoms is determined by only one (respectively two) scattering 
length difference(s). As can be seen from Tablelllll the param- 
eters Q,y and t5; nt extracted from a fit to the measured change 
of oscillation frequency versus magnetic field are directly re- 
lated to those scattering length differences. 

For / = 1 the difference (02 - «o) directly results from a fit 
of the form (|S) to the measured data. We perform this fit with 
one free parameter S\ , the interaction detuning for the process. 
The fit yields a value of 61 = 2k X (-11 .8^}° 4 ) Hz. We stress, 
that this parameter is still free from systematic effectes due to 
our uncertainty in optical potential depth, but has a statistical 
error, as well as a systematic error due to magnetic field uncer- 
tainty. The interaction strength Q,y can be expressed in terms 
of this parameter as Q,y = 2 y26\. 

For / = 2 the situation is more complicated. In partic- 
ular, the process |2, — 1;2, — 1) «-» |2, 0;2, -2) depends only 
on the scattering length difference (a\ - 02), whereas the dy- 
namics starting from |2, 0; 2, 0) depends on both differences 
(a a — a.2) and («2 - «o). We reanalyze the data presented in 
H33I1 . There, the dynamics was treated in a pure two-level 
model, and the Rabi-parameters were considered as indepen- 
dent. Here, we first fit the data set for oscillation frequency 
vs magnetic field of the process |2, — 1;2, — 1) «-> |2,0;2, -2) 
with only one parameter #2,2, reflecting the scattering length 
difference (a\ - 02)- The resulting best fit parameter is £2,2 = 
27r x (30.2*2-1) Hz. Then we fit the data set for the dynamics 
starting from |2, 0; 2, 0) with <52.2 fixed to this value, and leave 
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the second relevant interaction detuning 62,1 — reflecting the 
scattering length difference (02 - «o) — as fi" ee parameter. For 
this fit, which is shown in Figure^, we employ the analytical 
solution of the three level system Equations J13i . The best fit 
is achieved for 62,1 - 2nx (33.4*J'j) Hz. Due to the small un- 
certainty in each data point of the first data set compared to the 
uncertainty of the second data set, this treatment is equivalent 
to one combined fit of both data sets with two parameters. 

In order to extract the scattering length differences, we em- 
ploy Table [Hi] and the given values for U. Those scattering 
length differences are summarized in Table|l] The values now 
have an uncertainty due to the optical potential, entering via 
the constant U. 



f 


Aa 


Measured 


Predicted 


/=1 


(a 2 - a ) 


(-1.07^};{£ ±006 ) a B 


-1.38 a B 


/ = 2 


(a 2 - tfo) 


( 3 - 51 +(U8±0.16J fl B 


3.35 a B 


(fl4 — CI2) 


( 6 - 95 ±a07±0.26) fl B 


7.70 a B 



TABLE I: Scattering length differences inferred from the spin- 
oscillations measurement. The errors in the upper row give the sta- 
tistical uncertainty from the fit, whereas the first and second value in 
the lower row give the systematic error due to the uncertainty in mag- 
netic field and optical potential depth, respectively. The theoretical 
predictions are calculations based on values from 1381 14111 . 



D. Damping mechanisms 

An experimentally observed feature of the coherent spin dy- 
namics in both hyperfine states is the damping which is not ex- 
pected from the models describing the coherent dynamics. For 
the / = 2 manifold and high lattice depths the damping rate 
coincides with the loss rate of atoms due to hyperfine chang- 
ing collisions H33I1 . Although this loss process is absent for 
/ = 1, a finite damping on the order of y w 3 s -1 is observed 
even in the spin-1 data. One source for this damping could be 
the inhomogeneity of our system. Due to the gaussian inten- 
sity profile of the lattice laser beams, the wavefunction overlap 
J d 3 r\cf>o\ 4 is slightly higher in the center than in the outer re- 
gions of the system. Therefore the coupling strength Q,y and 
the constant detuning i5; nt are position dependent. This leads 
to a small dephasing in the system, with a rate that we esti- 
mate to be on the order of 0.3 s -1 . Therefore we can exclude 
its influence on the observed damping. Another possible ex- 
planation for the observed damping at high lattice depths in 
the spin-1 case are off-resonant Raman transitions introduced 
by the lattice laser beams. The expected rate of those events 
has been estimated to be similar to the observed damping rate. 
However, this effect has not been experimentally investigated. 
Also, heating to higher vibrational bands and subsequent tun- 
neling to neighboring sites is a possible damping mechanism. 

Another feature of the coherent spin dynamics in / = 1 
that has also been observed in f — 2 is an initial rise of the 
population in the final two-particle state, which is faster than 
predicted by the simple Rabi-model. This fast initial rise at 



the beginning of the spin dynamics together with the offset 
of the oscillations can partially be explained by the magnetic 
field ramp after the last microwave pulse. Here, the sample 
experiences a changing detuning, and the state from which 
the spin dynamics start from can be changed. A calculation 
including a model of the ramped magnetic field yields an off- 
set of the spin population oscillations which is a few 10% of 
the experimentally observed offset. 



V. NATURE OF THE MAGNETIC GROUND STATE 

As explained in the introduction, the relative orientation of 
two spins in a collisional event changes the actual interaction 
strength. In order to determine the particular orientation that 
minimizes the mean field energy at zero magnetic field, i.e. the 
magnetic ground state, it is convenient to express the interac- 
tion potential in terms of spin-independent and spin dependent 
parts (see e. g. IU2lll9ln . 

For the spin-1 case this yields the potential 



V s - co + c 2 /i ■ fi, 



(15) 



where cq = (4nh 2 /M) x (00 + 2«2)/3 is the spin-independent 
part, and C2 = (4nh 2 /M) x (02 - «o)/3 describes spin-spin 
interactions. For the antiferromagnetic or polar phase, the 
mean field energy according to Q is minimized by aligning 
the spins of two interacting atoms anti-parallel, and it emerges 
for C2 > 0, i.e. ai > ao 1131 . whereas the mean field energy 
in the ferromagnetic case is minimized by aligning the spins 
parallel, implying c<i < 0. 

In a similar way, the interaction potential for the spin-2 case 
can be written as 1 19] 



V s = c + ci/i ■ h + 5c 2 Po- 



(16) 



Here, Vq is a projector onto the singlet subspace, cq = 
(Anti 2 IM) x (4«2 + 3«4)/7 describes the spin-independent part, 
c\ = (4nh 2 /M) x (04 - «2)/7 determines the spin-spin interac- 
tion and C2 = (4nh 2 /M) x (7«o - 10«2 + 3«4)/7 accounts for 
the interaction between spin-singlet pairs. Differently from 
the spin-1 case, there exist three possible phases for a spin- 
2 Bose-gas at zero magnetic field. In addition to the ferro- 
magnetic and antiferromagnetic phases the system can be in a 
cyclic phase. In this phase the spin orientation of three inter- 
acting particles tend to form an equilateral triangle in a plane. 
The phase diagram is now two dimensional, depending on the 
spin-dependent coefficients C[ and ci- A sketch of the phase 
diagram is given in Figure |5] The ground state is ferromag- 
netic for c\ < and c\ — C2/2O < 0, antiferromagnetic for 
C2 < and c\ - C2/2O > 0, and cyclic for c\ > and ci > 0. 

A summary of the values for the spin-dependent coefficients 
in / = 1 and / = 2 determined from the frequency measure- 
ment of the coherent spin dynamics is given in Table |H] to- 
gether with the theoretically predicted values. For the spin-1 
case the inferred coefficient ci clearly shows a ferromagnetic 
ground state at zero magnetic field. For the spin-2 case the 
situation is not so clear. The coefficient c\ clearly excludes 
a ferromagnetic ground state, and the coefficient c^ indicates 
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FIG. 5: Phase diagram of the / = 2 magnetic ground state in 87 Rb. 



/ 


Ci 


Measured Aa/as 


Theoretical Ao/ob 


/ = 1 


C 2 /(47Th Z /M) 


-0 36 ±I,UI 


-0.46 


/ = 2 


cil(AnTi z IM) 


+0 99 ±U01 


+ 1.1 


c 2 /(4nh l /M) 


-0 53 ±(u 


-0.05 



tion[y]rernain unaffected by an overall change of U. The rea- 
son is that all arguments concerning the sign of the coefficients 
c\ and C2 hold equally for the measured interaction detunings 
<5i n t of the different processes, as long as the factor U is inde- 
pendent of the specific spin wave function. 

A second origin of the observed discrepancy between the 
experiment and the theoretical prediction might be that the 
description of the dynamics by only the spin-changing in- 
teraction potential Equation Q is not valid. For exam- 
ple, dipole-dipole interaction, although weaker than the spin- 
changing interaction, is also present in principle. In a per- 
fectly isotropic trap, they do not change the interaction energy. 
In an anisotropic lattice well, however, effects due to dipole- 
dipole interactions might affect the observed population os- 
cillations 1 54]. In this case, the conclusions drawn concern- 
ing the magnetic ground state of the two hyperfine manifolds 
might change H55I1 . This unresolved problem therefore needs 
further experimental and theoretical investigation. 



TABLE II: Summary of the measured coefficients of the spin- 
dependent interaction <15I16> . The errors are the statistical errors 
(upper row) and the systematic errors (lower row, first number due to 
magnetic field, second due to lattice potential). The calculated values 
are based on recently predicted scattering lengths 13a. 14111 . 



an antiferromagnetic ground state, but the error bars include 
the phase boundary to the cyclic phase. Moreover, we have 
observed that even slight changes in, e.g., the fit procedure 
can change the sign of cj- We conclude that although the data 
measured data strongly points towards an antiferromagnetic 
ground state, our analysis does not allow a final determination 
of the ground state in / = 2. 



VI. CONCLUSION 

In summary, we have shown that both hyperfine states of 
87 Rb show collisionally driven spin oscillations between two- 
particle states in a deep optical lattice. The oscillations can be 
described by a Rabi-model for a broad range of parameters. 
In the case where two final states are accessible, the observed 
spin dynamics can be explained by a coupled three-level sys- 
tem. 

The parameters of the Rabi model are directly related to dif- 
ferences of the scattering lengths in the 87 Rb hyperfine states 
/ = 1 and / = 2. The measured values can be compared 
to calculations based on recent theoretical predictions of the 
scattering lengths, and show good agreement. We observe, 
however, that a discrepancy on the order of 10% persists in all 
cases (see also 1 53]), which cannot be explained by either the 
experimental error bars or the uncertainty in the theoretically 
predicted scattering lengths. This discrepancy may have two 
reasons. First, our approach of using non-interacting wave 
functions as introduced in Section Hi Al in order to describe 
the wave function overlap might break down for deep lattices. 
This problem still requires further theoretical investigation. It 
should be noted, however, that the conclusions drawn in Sec- 



Acknowledgments 

We would like to thank Servaas Kokkelmanns and Johan 
Mentink for helpful discussions and for providing us with the 
s-wave scattering lengths for / = 2. This work was supported 
by the DFG, the EU (OLAQUI) and the AFOSR. FG acknowl- 
edges support from a Marie-Curie fellowship. 



COUPLING STRENGTHS 

It is important to notice that for the collision of two spin- 
1 atoms only one scattering length difference is sufficient to 
describe the dynamics; for the case of two colliding spin-2 
particles there are two relevant scattering length differences. 
For / = 1 this is the difference at - ao, and for / = 2 there 
are two differences at - ao and a\ - at- In most of the cases 
it is, however, more convenient to use derived quantities that 
directly describe the dynamics for the actual experimental sit- 
uation. 



Rabi dynamics 

Following the line of argument given in Section |H] the pa- 
rameters Q,y and 5; nt of the Rabi model can be calculated. Ta- 
ble|ni]summarizes the dependency of each relevant parameter 
in the Rabi model on the specific scattering length. It should 
be noted that there the scattering lengths ao and «2 have dif- 
ferent values depending on the hyperfine manifold, / = 1 or 
/ = 2. 

As described in Section HVl we fit the measured oscillation 
frequency vs magnetic field with only one (two) free param- 
eters for / = 1 (/ = 2), independent of our uncertainty in 
optical potential depth. For this, we express all relevant quan- 
tities in terms of interaction detunings 5; nt (see last column in 
TableED. 
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/ 


Process 


Parameter 


Corresponding Aa 




/ = 1 


|0, 0)^ | + 1,-1> 


Oi f 


2 V2/3(a 2 - «o) 


2V2.S, 


(5 int = S\ 


(a 2 - a )/3 


5, 


/ = 2 


|0, 0> <-> |+ 1,-1) 


a if 


2 V2/35[7(a 2 - a ) + 12(a 4 - a 2 )] 


2V2(<5 2 , 1 +2<5 2-2 ) 


Smt = S 2 ,l 


[l{a 2 - a ) + 2(a 4 - a 2 )]/35 


<?2,1 


/ = 2 


|0,0><-> | + 2,-2> 


Q., f 


2 V2/35[-7(a 2 - a ) + 3(a 4 - « 2 )] 


2V2(<5 2 , 2 -<5 2 , 1 ) 


6mt 


[7(a 2 - cto) + 17(a 4 - « 2 )]/35 


(5 2 ,i + 3(5 2-2 ) 


/ = 2 


1 + 1,-1)^1 + 2,-2) 


Ckf 


4/35[7(a 2 - «o) + 2(a 4 - a 2 )] 


4(5 2 , 2 


Smt 


3/7(a 4 - a 2 ) 


3(5 2 , 2 


/ = 2 


1-1,-1)^10,-2) 


a, f 


4V3/7(a 4 -a 2 ) 


4V3(i 2 , 2 


Smt = #2,2 


(« 4 - "2)/7 


<^2,2 



TABLE III: Summary of the collisional coupling strengths and detunings for the experimentally investigated cases as a function of the spin- 
dependent scattering lengths a F . The last column gives expressions of the various parameters in terms of interaction detunings (5, of the 
measured processes used in the fitting procedure predented in Section llV CI 
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